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Abstract 

We present a high precision Monte Carlo study of various universal am¬ 
plitude ratios of the three dimensional Ising spin model. Using state of the 
art simulation techniques we studied the model close to criticality in both 
phases. Great care was taken to control systematic errors due to finite 
size effects and correction to scaling terms. We obtain C+/C- = 4.75(3), 
f+, 2 nd/f-, 2 nd = 1-95(2) and u* = 14.3(1). Our results are compatible with 
those obtained by field theoretic methods applied to the </> 4 theory and high 
and low temperature series expansions of the Ising model. The mismatch 
with a previous Montecarlo study by Ruge et al. remains to be understood. 
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1 Introduction 


In the neighbourhood of a second order phase transition various quantities display 
a singular behaviour. In this limit most of the microscopic features which charac¬ 
terize a given model become irrelevant and models which differ at the microscopic 
level may share the same singular behaviour. This is the basis of the concept of 
universality. The first, well known, consequence of universality is that different 
models belonging to the same universality class share the same critical indices. 
However the hypothesis of universality of the various scaling functions has much 
stronger implications and it is possible to show that models belonging to the same 
universality class are also characterized by the same values of some critical-point 
amplitude combinations |I|. Let us see a simple example. Near the critical point 
the correlation length £ diverges as 

£ ~ 

£ ~ 

with 


where T is the temperature and T c the critical temperature. Different models in 
the same universality class share not only the same critical exponent u, but also the 
same dimensionless combination of critical amplitudes /+//_. This is particularly 
relevant from the experimental point of view, since in general critical amplitudes 
are more easily detectable than critical indices and allow a simpler identification of 
the universality class. In fact the variations of the critical indices between different 
universality classes are in general rather small, while the amplitude ratios may 
vary by large amounts. In this paper we shall be in particular interested in the 
universality class of the three dimensional Ising model which has several interesting 
experimental realizations, ranging from the binary mixtures to the liquid vapor 
transitions. 

The two standard approaches to the evaluation of these amplitudes ratios in the 
Ising case are the use of field theoretic methods applied to the £> 4 theory 0-||, 
and the extrapolation to criticality of low and high temperature series expansions 
on various lattices |D|]. All these estimates are in general in rather good agreement 
among them (for a comparison and a discussion see sect.5). 

In order to obtain results from Montecarlo simulations relevant for the scaling 
limit we have to control both finite size effects as well as corrections to scaling. 
This means that the linear lattice sizes L have to be chosen such that L » £ while 
£ —> oo as f3 —> /3 C . In practice one has to carefully check which factor of L/£ is 
required to obtain results sufficiently close to the thermodynamic limit. While in the 
high temperature phase L/C, ^ 7 turns out to be sufficient to give thermodynamic 
limit results within numerical accuracy, in the low temperature phase this factor 
has to be doubled at least. The value of £ that can be reached, and hence the 


f+ t t > 0 
t< o 


t = 


T-T r 


( 1 ) 

( 2 ) 
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control of corrections to scaling, is limited by the CPU time available for the study. 
In the present paper the largest correlation length is £ = 11.884(9) in the high 
temperature phase and £ = 6.208(18) in the low temperature phase. 

One also should note that simulations in the low temperature phase are consid¬ 
erably more difficult than those in the high temperature phase of the model. In the 
low temperature phase conceptual as well as practical problems caused by sponta¬ 
neous symmetry breaking arise. Furthermore the determination of the correlation 
length is complicated by the occurrence of secondary correlation lengths which are 
close to the leading one. 

For all these reasons only in these last years there have been some attempts 
to measure these ratios in Montecarlo simulations 0. However the results are 
rather puzzling. For instance, in the case of the ratio /+//_ discussed above, the 
Montecarlo estimate, which is / + //_ = 2.06(1) ]nj], disagrees with the one obtained 
with strong/weak coupling series /+//- = 1.96(1) |9], while the field theoretical 
estimate /+//_ = 2.013(28) || lies in between the two. 

The aim of our work is to show that Montecarlo estimates of the amplitude ratios 
can indeed be competitive with other approaches. To this end we have elaborated 
a technique to directly extract the various amplitude ratios, without evaluating the 
single amplitudes thus avoiding all the uncertainties related to the critical indices. 
We shall discuss this point in sect.5 below. Besides this, we have devoted a great 
care throughout the paper to keep under control systematic errors due to finite 
size effects and corrections to scaling. Finally we have used state of art simulation 
techniques to obtain high precision estimates of the observables near the critical 
point, in both phases. The simulations in the high temperature phase have been 
performed using Wolff’s single cluster algorithm. Here the improved estimators give 
a great boost to the accuracy of the results. However in the low temperature phase, 
due to the finite magnetization, the improved estimators of the cluster-algorithm 
are of little help. Hence we simulated here with a multispin coding implemented 
Metropolis-like algorithm. 

As we shall see our results are comparable in precision and agree with the most 
recent field-theoretic and strong/weak coupling estimates, while they are incompat¬ 
ible with the MC results of Ruge et al. 0 - A more detailed comparison of the 
data might be helpful to understand this discrepancy. 

This paper is organized as follows: In sect.2 we have collected some informations 
on the three dimensional Ising model and on the observables that we study. In 
sect.3 we discuss the details of the simulation, while in sect.4 we analyse the scaling 
behaviour of the measured quantities: magnetization, susceptibility and correlation 
lengths. In sect.5 we study the amplitude ratios and compare our results with 
other existing estimates and with the experiments. Finally sect.6 is devoted to 
some concluding remark. 
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2 General Setting 


2.1 The Model 


We study the Ising spin model in three dimensions on a simple cubic lattice. The 
action is given by 


Sspin V 


E 


Sn^r, 


( 3 ) 


<n,m> 


where the held variable s n takes the values —1 and +1; n = (no,rii,n 2 ) labels 
the sites of the lattice and the notation < n, m > indicates that the sum is taken 
on nearest neighbour sites only. The coupling f3 is defined as f3 = hence the 
reduced temperature t can be written as 


Pc~ P 


( 4 ) 


where P c = -r^r. We shall consider in the following ri\ and n 2 as “space” directions 
and no as the “time” direction and shall sometimes denote the time coordinate Uq 
with r. We always consider lattices of equal extension L and periodic boundary 
conditions in all three directions. 


2.2 The observables 


2.2.1 Magnetisation 

The magnetization of a given configuration is defined as: 


771 


v^ Si 1 


( 5 ) 


where V = L 3 is the volume of the lattice. However, in a finite volume the Z 2 
symmetry of the model can not be broken for any nonzero temperature. Hence the 
expectation value of m vanishes. 

In order to obtain the magnetization of the model in the low temperature phase 
one should add a magnetic held h in order to break the symmetry. Then one should 
first take the thermodynamic limit at finite magnetic held and then take the limit of 
vanishing magnetic held. However it is difficult to follow this route in a numerical 
study. 

As an alternative Binder and Rauch 


12| suggested to simulate the finite lattices 


at vanishing external held and study the quantity 


< m > = lim V< 77i 2 > 


( 6 ) 


However it turns out that this is not the best choice. In fact this observable is 


affected by strong finite size effects [13 which would require very large lattices to 
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obtain reliable estimates of the infinite volume magnetization. It has been recently 


observed [14] that a much more stable observable is: 


< m > 


= lim < 

L —>oo 


m\ > 


( 7 ) 


The finite size behaviour of this observable, was carefully studied in [|14|] where 
it was shown that the asymptotic, infinite volume, value is reached for lattices of 
size L >~ 8£, where £ denotes the correlation length. In our simulations we always 
used lattice sizes much larger than this threshold. 

Close to the critical temperature, the magnetization is supposed to scale as 


< m >~ B (— t) 13 ] t < 0 


( 8 ) 


where the critical exponent [3 should not be confused with the inverse temper¬ 
ature. 


2.2.2 Magnetic susceptibility 

The susceptibility gives the response of the magnetization to an external magnetic 
held. 


d < m > 

dli 


( 9 ) 


One easily derives that the magnetic susceptibility can be expressed in terms of 
moments of the magnetization as follows: 

X = v(<m 2 > — <m > 2 ) . (10) 

Close to the critical temperature the magnetic susceptibility is supposed to scale as 


X ~ C+ t 7 ; t > 0 

x ~ CCH)- 7 ; t< o . (li) 


2.2.3 Exponential correlation length. 

We consider the decay of so called time-slice correlation functions. The magnetiza¬ 
tion of a time slice is given by 


Sn 0 j 2 S (n 0 ,ni,n 2 ) 

^ 711,712 


( 12 ) 


Let us define the correlation function 

GW=y{(S„ 0 S„„ +T )-(S„ 0 ) 2 } 

n 0 

The large distance behaviour of G(r) is given by 

G(t) oc exp(— t/£) , 


(13) 


(14) 
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where £ is the exponential correlation length. 

Close to criticality the behaviour of the correlation length is governed by the 
scaling laws 


£ ~ /+ t u ; t > 0 

£ ~ /_ H)H t < o 

(15) 

2.2.4 Second moment correlation length. 


The square of the second moment correlation length is defined for 
of the spacetime dimensions d by 

a generic value 

C 2 /-D 

?2nd 2d/i 0 ’ 

(16) 

where 


/Jo = hm — V (. s m s n ) c 

L—> oo V 

v m,n 

(17) 

and 


h 2 = hm — V (m - n) 2 (s m s n ) c . 

L— XX) 1/ 

v m,n 

(18) 

The connected part of the correlation function is given by 


(SmSn)c (km^n) (^m) 

(19) 


This estimator for the correlation length is very popular since its numerical 
evaluation (say in Montecarlo simulations) is simpler than that of the exponential 
correlation length. Moreover it is the length scale which is directly observed in 
scattering experiments. However it is important to stress that it is not exactly 
equivalent to the exponential correlation length. The relation between the two can 
be obtained as follows. Let us write 

/i 2 = E ( n ~ 171 ) 2 ( SmS n )c 

v m;n 

^ d— 1 

= 77 E E K “ "V) 2 i S mSn)c 

* n;m ft =0 

TT ) ) (no m o) (SmSn)c ■ (20) 

' n;m 

Due to the exponential decay of the correlation function this sum is certainly con¬ 
vergent and we can commute the spatial summation with the summation over con¬ 
figurations so as to obtain 

OO 

^2 = d E ^ (So S T ) C (21) 
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with S no given by eq. ([T^) . Analogously one obtains 


OO 

hO = E (^O^r)c 

T— — OO 


( 22 ) 


If we now insert these results in eq.([L 6 ]), assume a multiple exponential decay 


(S 0 S r ) c ocJ2 c i exp(-|r|/&) , 


and replace the summation by an integration over r we get 

,2 1 /r”o dTT2 expl-T/O Ei Cjff 

i2nd 2 /“ 0 dr exp(-T/0 E,c,5, 


(23) 


(24) 


which is equal to £ 1 2 if only one state contributes. An interesting consequence of 
this analysis is that the difference from one of the ratio £,/&nd gives an idea of the 
density of the lowest states of the spectrum. If these are well separated the ratio 
will be almost one, while a ratio significantly higher than one will indicate a denser 
distribution of states. 

The critical behaviour of £2 nd is governed by the same critical index v, so near 
the critical point we expect: 


^ 2 nd ~ f+,2nd t , t 0 

6 nd ~ f-,2nd(-t)- U ] t< 0 (25) 


2.3 Critical indices 


Our aim is to obtain high precision estimates for some amplitude ratios. To this end 
we need to use as input informations the critical temperature (3 C and the values of 
the critical indices Q defined above. We list in tab.l some estimates for these quanti¬ 
ties, obtained with held theoretical methods, strong coupling series and Montecarlo 
simulations. 

In general these values are in rather good agreement among them, despite the 
fact that they were obtained with very different methods. As input parameters for 
our analysis we have decided to choose the following values: 


P c = 0.2216544(3) , 7 = 1.2390(15) , v = 0.6310(15) , (26) 

P = 0.3270(6) , 6 = 0.51(3). (27) 

which are obtained combining together the results of |b|] and [|l 6 j and satisfy the 
scaling relations. Let us stress however that our results are only slightly affected 
by this choice |j|. 


1 As a matter of fact, for the actual determination of the amplitude ratios we only need to know 
/3 C and 9. The values of the other critical indices will only be used in comparing our results with 
those obtained with the series expansions. 

2 The systematic errors due to the uncertainties in the choice of (3 C and 9 turn out to be much 
smaller than the statistical fluctuations of our estimates. In any case, both are taken into account 
in the errors that we quote in our final results. 
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Table 1: Results for /3 C and for the critical indices given in the literature 


ref. 

method 

Pc 

7 

V 

P 

6 


R 


FT 


1.241(4) 

0.630(2) 

0.324(6) 

0.496(3) 


16 


e-expansion 


1.2390(25) 

0.6310(15) 

0.3270(15) 

0.51(3) 


17 


d = 3 


1.2405(15) 

0.6300(15) 

0.3250(15) 

0.50(2) 


18 


MCRG 

0.221652(3) 


0.624(1) 


0.50-0.53 


19 


MCRG 

0.221655(1)(1) 


0.625(1) 


0.44 


2 (J 


MC,FS 

0.2216546(10) 

1.237(2) 

0.6301(8) 

0.3267(10) 

0.52(4) 


21 


MC 

0.2216576(22) 

1.239(7) 

0.6289(8) 

0.3258(44) 



14 


MC 

0.2216544(3) 



0.3269(6) 

0.508(25) 


22 


HT 


1.237(2) 

0.6300(15) 


0.52(3) 


23 


HT 


1.239(3) 

0.632(3) 


0.52(3) 


24 


HT 


1.2385(25) 

0.6305(15) 


0.57(7) 


25 


HT 


1.2395(4) 

0.632(1) 


0.54(5) 


2.4 Amplitude ratios 

In the following we shall be interested in these scaling functions: 


r x (f) 


x(t) 
x(-t) ’ 


r &) 


f,2 nd,{t) 
£2 nd(~t) 


(t > 0 ), 


u{t) 

r e (t) 


3 x(t) 
tLdit)™ 2 ® 
x(t) 




(t < 0 ), 

(t > 0 ) 


(28) 

(29) 

(30) 


(notice the factor of three difference between the definitions of T c and u). While 
T x , T? and r c mix low and high temperature observables, u only contains quantities 
evaluated in the broken symmetry phase. T c and u are scale invariant thanks to 
the following scaling (and hyperscaling) relations among the critical exponents: 


ct + 2/3 + 7 = 2 , dis = 2 — a 


(31) 


In particular, u plays the important role of a low temperature renormalized coupling 
constant in the study of the </> 4 theory directly in d = 3. 

It is important to notice that T c is related to the ratio of two amplitude combi¬ 
nations (in which A(t) denotes the specific heat): 

Rc = Xit) v A( y R^bndWAit) 1 ' 3 (t> 0) (32) 

which have been widely studied in the literature since they can be rather easily 
evaluated in experiments. The relation is: T c = R c /R |. In the scaling limit these 
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functions are related to the amplitudes defined in eq. 


25|) as follows: 


fes = it • !“S r «(‘) = yy) 


!'h “<*> = = 75^ • !‘h r <« = 


Ch 


f 3 R 2 ' 
J + ,2nd n 

Finally we shall also be interested in evaluating the ratio: 

£ 


£2 nd 


(33) 

(34) 

(35) 


both above and below the critical point. 


2.5 Series expansions 

A very powerful approach to the study of the three dimensional Ising model is 
represented by the series expansions which lead to estimates for several quantities 
near the critical point which are competitive with the most precise Montecarlo 
simulations. I 11 the following we shall compare our results in the low temperature 
phase (which is the one in which simulations are more difficult and results are 
in general affected by stronger finite size effects) with those obtained with series 
expansions with the twofold aim of testing the reliability of our simulations and 
of comparing the precision of the two methods. The low temperature regime is 
also particularly interesting because recently these series have been extended up 
to very high orders |26], ^7], |28|. Some informations on these series can be found 


in tab. 13. In order to extract from the series the estimates for the observables in 
which we are interested and to quantify the uncertainty of such estimates we use 
the so called “double biased inhomogeneous differential approximants” (IDA). The 
technique of IDA is described in ||, |29|, to which we refer for notations and further 
details. Following ref. J9|] we use the notation [K/L;M] for the approximants. In 
order to keep the fluctuations of the results under control, we have chosen to use 
double biased IDA namely we fix both the critical coupling /3 C and the critical 
index describing the critical behaviour of the observable. As K and L vary we 
obtain several different IDA’s and correspondingly several different estimates of the 
observable. As we approach the critical point these estimates start to spread out, 
indicating that we are pushing the series toward its convergence threshold. The 
last problem is then to extract from this set of values the best estimate and its 
uncertainty. Our choice in this respect is to neglect those IDA’s which fluctuate 
too wildly and treat the remaining approximants on the same ground. To this end 
we determine the smallest interval that contains half of the results. The values 
that we shall quote in the tables below as our best estimates correspond to the 
centre of this interval, the first number in brackets gives the half of the size of the 
interval. The second number in brackets gives the error induced by the error of 
















(3 C = 0.2216544(3) and of the critical index used for biasing. Together they give an 
idea of the uncertainty of the estimate. As we shall see, in the range of coupling in 
which we are interested, the uncertainty will be always dominated by the spread of 

IDA’s. 


3 The simulations 


3.1 Simulations in the low temperature phase. 

We simulated the Ising spin model in its low temperature phase at f3 = 0.2391, 
0.23142, 0.2275, 0.2260,0.2240 and 0.22311 using a demon-algorithm implemented 
in the multispin coding technique. A detailed discussion of this algorithm is given 
in ref. 0. The update of a single spin takes 46 x 10 9 sec. on a HP 735 and 


21 x 10 ~ 9 sec. on a DEC Alpha 250 workstation for L = 120. For (3 = 0.22311 the 
integrated autocorrelation time of the magnetization was Ti nt = 81.(2.) in units of 
sweeps. 

We used cubical lattices with periodic boundary conditions and a linear exten¬ 
sion of about 20£. It should be noted that test runs revealed that in contrast to the 
high temperature phase a linear lattice size of 6 £ is clearly not sufficient to obtain 
results close to the thermodynamic limit. Some informations on the simulations are 
collected in tab. 2 . 

We computed {So S T ) for all values of r available on the finite lattice. We 
evaluated Sq S t for all three lattice directions and all possible translations. The 
connected part was then obtained by subtracting the expectation value of the square 
of the magnetization {{1/V X)„s n ) 2 ). 

In order to obtain an estimate for the true correlation length we started from 
the ansatz 

G(t) oc exp(— t/£) + exp(—(L — r)/£) , (36) 

where the last term takes into account the periodicity of the lattice. An effective 
correlation length £ e //('r) is then computed by solving the equation above for r and 
r + 1. Ignoring the term exp(— (L — r)/£) , £ e ff( r ) takes the form 

^ e//(r) = ln(G(r + l)) -ln(G(r )) 5 ^ 

For r > 3£ £ e //( r ) seems to stabilize within error bars. In tab.3 the results for 
r = 3£ are given. However it is important to stress that there might still be 
systematic errors due to higher excitations that are of the same magnitude as the 
statistical error of £ e //- Therefore a multi-exponential ansatz might be useful. We 
shall further discuss this point in sect. 4.4 below. 

We computed the second moment of the correlation function by 

Tmax OO 

/A2 = 51 t2G \ t ) + t2 G '( T max) exp( — (t - T max ) / £ e ffurnace)) (38) 

T= 1 'T=Tmax 1 
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where G'{r ) = C e ff{r) exp(—r/^ e //(r)). Where again C e ff{r ) and £ e //(r) are 
obtained from 

G(t) oc exp(-r/0 + exp(-(L - r)/f) , (39) 

inserting r and r + 1. As for the exponential correlation length we used r max = 3£ 
for the data reported in tab.3 . Note that the systematic error introduced by the 
finite T max only affects the second term of eq. |3^, which is small compared to the 
first one. Hence these systematic error can safely be ignored for our choice of T max . 
The susceptibility is computed analogously. The results are summarized in tab. 3. 


Table 2: Statistics of the runs in the low temperature phase. 


p 

L 

measures 

sweeps/measure 

bits 

0.2391 

30 

40000 

25 

64 

0.23142 

40 

50000 

25 

64 

0.2275 

50 

50000 

25 

64 

0.2260 

80 

50000 

25 

64 

0.2240 

100 

92000 

25 

32 

0.22311 

120 

124000 

25 

32 


Table 3: Results in the low temperature phase 


P 

L 

m 

E 

S,exp 

£,2nd 

X 

0.2391 

30 

0.667162(20) 

0.553732(17) 

1.2851(28) 

1.2335(15) 

4.178(3) 

0.23142 

40 

0.570306(16) 

0.478046(12) 

1.8637(45) 

1.8045(21) 

9.394(4) 

0.2275 

50 

0.491676(14) 

0.430364(10) 

2.578(7) 

2.5114(31) 

18.706(10) 

0.2260 

80 

0.449984(16) 

0.409609(4) 

3.103(7) 

3.0340(32) 

27.596(11) 

0.2240 

100 

0.372490(10) 

0.378612(3) 

4.606(13) 

4.509(6) 

61.348(34) 

0.22311 

120 

0.320830(10) 

0.362946(2) 

6.208(18) 

6.093(9) 

112.60(7) 


3.2 Simulations in the high temperature phase. 


We simulated the Ising spin model in the high temperature phase using the single 
cluster algorithm |3T| . The time slice correlation function was determined using the 
cluster improved estimator. Finite size effects are less important in this phase and 
some preliminary test showed that lattice sizes greater than 6£ are enough to keep 
them under control. This is clearly visible in the data of tab.4 where we report a 
test at (3 — 0.21931. 

Also the determination of the correlation length in the high temperature phase 
turns out to be much easier than in the low temperature phase. The effective cor¬ 
relation length approaches a plateau quite quickly and the true correlation length 
is well approximated by at a self-consistently chosen distance r ~ £ e ff. This 
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Table 4: Test for finite size corrections at (3 = 0.21931. 


p 

L 

stat 

f,exp 

£,2nd 

E 

X 

0.21931 

40 

50000x1000 

8.701(4) 

8.691(5) 

0.313986(24) 

303.23(30) 

0.21931 

50 

50000x1000 

8.747(4) 

8.741(5) 

0.313870(18) 

307.14(31) 

0.21931 

60 

50000x1000 

8.754(6) 

8.750(5) 

0.313811(15) 

307.79(31) 

0.21931 

70 

50000x1000 

8.758(4) 

8.751(5) 

0.313823(14) 

307.95(31) 

0.21931 

80 

50000x1000 

8.766(5) 

8.760(5) 

0.313849(14) 

308.58(31) 


behaviour of £ e // also implies that the difference between the second moment cor¬ 
relation length and the true correlation length is much smaller that in the low 
temperature phase. We have chosen the (3 values in the high temperature phase 
such that (3 C — f3 = /3i ow — (3 C , where /3i ow are the inverse temperature used in the 
simulations of low temperature phase. The reason for this choice will be clear in 
the following section. The uncertainty of (3 C virtually does not affect the following 
analysis. The simulation results are summarized in tab.5. Recently in ref. 


Table 5: Results in the high temperature phase 


P 

L 

stat 

£exp 

S2 nd 

E 

X 

0.20421 

20 

50000x500 

2.363(1) 

2.346(1) 

0.262928(49) 

25.255(16) 

0.21189 

30 

50000x500 

3.477(2) 

3.465(2) 

0.284663(35) 

52.09(4) 

0.21581 

40 

50000x500 

4.864(3) 

4.854(3) 

0.298366(28) 

98.90(10) 

0.21731 

50 

50000x800 

5.892(3) 

5.885(3) 

0.304493(20) 

143.04(12) 

0.21931 

80 

50000x1000 

8.766(5) 

8.760(5) 

0.313849(14) 

308.58(31) 

0.22020 

100 

50000x1200 

11.884(9) 

11.877(7) 

0.318742(11) 

557.57(61) 


Monte Carlo results for the second moment correlation length and the magnetic 
susceptibility were reported. Interpolation of their results, using the scaling ansatz, 
to our /^-values leads to results consistent with ours. One has to note however, 
that our statistical errors are at least 3 times smaller than those of ref. [[hj in the 
common /3-range. 


4 Analysis of the results 


4.1 Magnetization 

A very precise Montecarlo study of the behaviour of the magnetization in the Ising 


model can be found in |TJJ. In particular, in JT4J it was shown that the magneti¬ 


zation values obtained from the montecarlo simulations were well described by the 
following empirical approximation: 


m(l 3) = f 


0.32694109 


(1.6919045 - 0.34357731 f 


0.50842026 


- 0.42572366 t) 


(40) 
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with t given by eq.([4j). Since our values of /3 are inside the region of validity of 
this approximation, it is interesting to compare also our magnetization values with 
eq.(fltj). Notice as a side remark, that our estimates for the magnetization are in 
general more precise of those reported in [|b| (where however a much larger number 
of (3 values was studied). This comparison is reported in tab.6, together with the 


estimates obtained with a double biased IDA analysis of the series published in p8 


Table 6: Comparison of our Monte Carlo results for the magnetization with eq. |4C 
and with double biased IDA’s. 


(3 

our MC 

eq. 10 of I4| 

biased IDA 

0.2391 

0.66716(2) 

0.667143 

0.667151(3)(1) 

0.23142 

0.570306(16) 

0.570279 

0.570300(16)(1) 

0.2275 

0.491676(14) 

0.491645 

0.49167(4)(1) 

0.2260 

0.449984(16) 

0.449953 

0.44999(7)(1) 

0.2240 

0.372490(10) 

0.372471 

0.37253(15)(2) 

0.22311 

0.320830(10) 

0.320809 

0.3209(2)(1) 


It is interesting to notice that our values are always in perfect agreement with 
those obtained from the series expansions, and that our results become more precise 
than the strong coupling ones starting from £ ~ 2. 


The agreement with the data of ref jl4| is also very good. Even if our values are 
systematically slightly higher than those of eq. (|40|) , they are well inside the error 
bars reported in tab.l of ref |14| . 


4.2 Susceptibility 

In tab.7 we report the comparison of our data on the susceptibility in the low 
temperature phase with a double biased IDA analysis of the series published in |ES 


The agreement is again very good and as in the previous case our results become 
more precise than the strong coupling ones starting from £ ~ 2. 


Table 7: Comparison of our Monte Carlo results for the susceptibility with double 
biased IDA’s. 


/3 

our MC 

biased IDA 

0.2391 

0.23142 

0.2275 

0.2260 

0.2240 

0.22311 

4.178(3) 

9.394(4) 

18.706(10) 

27.596(11) 

61.348(34) 

112.60(7) 

4.1801(16) 

9.401(20) 

18.76(15) 

27.67(40) 

61.8(2.7) 

114.6(10.5) 
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4.3 Second moment correlation length 


In tab.8 we report the comparison of onr data on the susceptibility in the low 


temperature phase with a double biased IDA analysis of the series published in |26 


Also in this case the agreement is very good. In addition we give Montecarlo results 
33j for which the /3-value matches with ours. One has to note that the results 

4.7£ and L 


of ref. 
of ref. 


|33| were obtained with L 
respectively. 


7.8£ for 0 = 0.226 and j3 = 0.224 


Table 8: Comparison of our Monte Carlo results for the second moment correlation 
length with double biased IDA’s. 


p 

our MC 

MC of ref 

0 

biased IDA 

0.2391 

0.23142 

0.2275 

0.2260 

0.2240 

0.22311 

1.2335(15) 

1.8045(21) 

2.5114(31) 

3.0340(32) 

4.509(6) 

6.093(9) 

3.22(1) 

4.61(6) 

1.2358(16) 

1.803(5) 

2.509(11) 

3.034(16) 

4.493(30) 

6.084(46) 


4.4 Exponential correlation length 


As we mentioned above, in the low temperature phase the evaluation of the expo¬ 
nential correlation length is much more delicate than in the high temperature phase. 
In particular, we know from the fact that the ratio 7-^— is significantly different from 

s2 nd 

1 and from [Q that in this region the spectrum is very rich and that nearby states 
exist that could contaminate the measure of £. This is exactly the situation dis¬ 
cussed in sect.3.1 above and accordingly we may expect some systematic error in 
ieff- Since the presence of nearby masses is a rather common situation in the bro¬ 
ken symmetry phases of statistical mechanical models and since, notwithstanding 
this, the estimator £ e yy is commonly used also in this cases, we have decided to 
devote this section to a detailed analysis of this problem. We can explicitly see that 
£ e ff evaluated according to eq.fl37|) is not a good estimator of the true correlation 
length by comparing our estimates with those of ref. [Q (see tab.9). In ref. |HJ we 
computed the glueball spectrum of the Z 2 gauge theory in 3 dimensions. In d = 3 
the spin and gauge Ising models are related by duality and the inverse of the 0 + 
glueball mass exactly coincides with exponential correlation length of the spin Ising 
model. In ref. [|34| we used a variational approach, using 27 different wilson loops 


as operators, to obtain a faster convergence of £ e ff- bi this way each mass of the 
spectrum was driven in a different channel, and practically no contamination from 
higher states was present. The results of ]R4| are comparable in statistical accuracy 
with the ones presented here. It is easy to see looking at tab.9 that the values of 
£ e ff obtained here are systematically smaller, in average by a factor of 0.995(2). 
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This shows, as expected, that the single exponential ansatz is problematic in this 
case and that a multi-mass ansatz is needed. In order to have an independent test 
of this fact we tried to fit our data at (3 = 0.224 for the correlation function with 
the following 3-mass ansatz 

G{t) ~ ci exp(— t/£l) + c 2 exp(— t/£ 2 ) + c 3 exp(-r/£ 3 ) . (41) 


However, the main problem of such multi-mass fits is that they are in general 
rather unstable under variation of the fit-range. This was also the case with our 
fits. Therefore we fixed the values of £ 2 = 2.50 and £ 3 = 1.70 found in ref. [34 . 
We found the following values: c 2 /ci ~ 0.1 and c 3 /ci ~ 0.04, when r’s in the range 
of 5 to 25 are included in the fit. Notice however that even in this case the results 
were still rather unstable, and for this reason we cannot give reliable error bars for 
our estimates. Assuming that G(t ) is well described by the three mass-ansatz and 
using our estimates for c 2 /ci and c 3 /ci we obtain £ e //(3£) = 0.993£ which is indeed 
consistent with our result above. 

It is also interesting to insert our estimate for c 2 /ci and c 3 /ci into eq. |24|. We 
obtain £/£ 2nd = 1-024. 

The situation is much simpler in the high temperature phase, where £ 2n d ~ £, 
no nearby masses are present and £ e ff is a good estimator of the true correlation 
length. 


Table 9: Comparison of the results for the exponential correlation length with those 
obtained for the Z 2 gauge model ((3 denotes the dual of (3) 


p 

P 

£gauge 

t*ff 

0.72484 

0.23910 

1.296(3) 

1.2851(28) 

0.74057 

0.23142 

1.864(5) 

1.8637(45) 

0.74883 

0.22750 

2.592(5) 

2.578(7) 

0.75202 

0.22600 

3.135(9) 

3.103(7) 

0.75632 

0.22400 

4.64(3) 

4.606(13) 


5 Universal amplitude ratios 

The standard approach to evaluate the amplitude ratios is to fit the data obtained 
for both phases separately with the expected scaling law, and then take the ratio 
of the amplitudes obtained from the fits. 

However the bias introduced by the uncertainty in the critical exponent can be 
avoided by directly studying the ratio as a function of the reduced temperature. 
This is the reason for which we carefully chose the couplings so as to have the same 
differences A/3 = \(3 — (3 C \ in the two phases. To better explain our approach let 
us study as an example the amplitude ratio T x . From the data reported in tab.3 
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and 5 we can compute the ratios of susceptibilities in the low and high temperature 
phase as a function of A/3. 


r x (A/3) 


x{Pc - A/3) 
xiPc + A/3) 


(42) 


As A/3 goes to 0 we expect T x (A/3) to converge to the amplitude ratio C + /C_. 
However the approach to this critical value is rather non-trivial. A naive implemen¬ 
tation of the scaling hypothesis would suggest that the data thus obtained should 
be constant within the error, but it is easy to see, by looking at the data in tab. 10 
that this is not the case. There are in fact two source of corrections. The fact that 
observables in both phases are involved in the ratio tells us that we must expect 
a correction proportional to A/3. Moreover we certainly expect a “correction to 
scaling” contribution proportional to A/3 e . In the example of T x the need of such 
corrections is clearly evident. Looking at the data in tab. 10 we see that the viola¬ 
tions of scaling are much larger than our statistical errors. Even for the smallest 
values of A/3 we see no stabilization of T x within error bars. Following the above 
discussion we fitted the data of tab. 10 with the law 


T x (A/3) = C+/C. + a 0 A [3 e + ai A/3 . (43) 

where we assumed that there is no other correction to scaling exponent 6' between 
9 and 1. The results of these fits for the various ratios in which we are interested are 
reported in tab.11. The fact that we always find rather low xf. ed strongly supports 
the correctness of the above assumption. 


Table 10: The various ratios as functions of A/3 


A/3 

r x 

r € 

u 

r c 

AT 

(.2nd 

0.01745 

6.044(5) 

1.902(2) 

15.00(6) 

4.394(9) 

1.042(3) 

0.00977 

5.546(5) 

1.920(3) 

14.75(5) 

3.850(10) 

1.033(3) 

0.00585 

5.283(3) 

1.932(3) 

14.66(6) 

3.577(10) 

1.026(4) 

0.00435 

5.182(5)(1) 

1.939(3) 

14.64(5) 

3.466(8) 

1.0227(34) 

0.00235 

5.027(6)(2) 

1.942(3) 

14.47(6) 

3.308(9) 

1.0215(42) 

0.00146 

4.947(6)(3) 

1.948(3)(1) 

14.51(7) 

3.233(9)(1) 

1.0188(45) 


Few comments are in order at this point. 

al] When dealing with combinations of observables all in the same phase we don’t 
need a correction to scaling term proportional to A/3. This is the case of the 
coupling constant u and of the ratio £/£2 nd- In this case we fitted with the 

lawf]: 

u( A/3) = u* + a 0 A/3 0 (44) 

3 Notice however that, hidden in the A/3 correction there should also be a term proportional to 
A/3 2S which, due to the fact that 9 ~ 1/2 is essentially indistinguishable from A/3. The correction 
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a2] The error due to the specific choice of /3 C is always very small. We have reported 
its value in the data of tab. 10 only when it is not negligible. In these cases the 
number in the first bracket gives the statistical errors of our data, while the 
second takes into account the uncertainty of the inverse critical temperature. 

a3] In the last row of tab. 11 we have reported onr final results for the various 
ratios. The corresponding errors take also into account the uncertainty in the 
index 6 and are thus slightly larger than those extracted by the fits. 


a4] In the results reported in tab. 11 we always fitted only the last five values of 
A/3, and systematically discarded the data at A/3 = 0.017450. 


a5] A particular care must be devoted to the study of the ratio . It is possible 
to prove that if higher masses exist in the theory (and this is the case in both 
phases of the Ising model) then must certainly be larger than 1. This is 
a consequence of eq.(^4|) and of the fact that the coefficients c t which appear 
in it must be positive (see eq. (9) and (10) of ref. [Q for a proof of this last 
statement). 

In the high temperature phase 7 -^— is almost compatible with 1 (within the 

s2 nd 

errors), and we can only use our data to set an upper bound for its value 
which, looking at the data with the largest correlation length, can be safely 
chosen to be f < 1.0006. 

J-\-,2nd 


On the contrary in the /3 > /3 C phase the quantity -1 is much larger than 

s2 nd 

the error bars, and can be measured rather precisely. The data in the last 
column of tab .10 and 11 refer to this case and use £ e // (defined in sect.4.4) as 
estimator of £. Hence we must add to the result of the fit (which is reported 
in the last row of tab.11: j, - = 1.009(5) ) the contribution due to the 
systematic underestimation A£ ~ 0.007 discussed in sect.4.4 above. Taking 
into account also this correction we quote as our final result = 1.017(7). 

It is interesting to notice that this result agrees within the errors with the 
value f~ d ~ 1.024 obtained in sect. 4.4 by inserting in eq.(^4|) our estimates 
for C 2 /C 1 and C 3 /C 1 and the values of the two nearby (inverse) masses £2 and 


A [3 2S should be present also in the case in which all the observables belong to the same phase, 
thus suggesting to use also in this case the fit eq.(fl3|). It turns out however that such a A (3 26 
correction, if present, has a negligible amplitude and for this reason we confined ourselves to the 
fit (@. 

4 This is not due to the fact that adding this data we had a poorer fit, on the contrary we 
checked that that for all the ratios the fit keeping all the six data was always equally good. The 
reason of our choice is that we tried to confine ourselves to the narrowest possible region near the 
critical point compatible with a reasonable precision for the results. This allows us to trust in 
our assumption of neglecting other possible, unknown, corrections to scaling which are certainly 
present but hopefully negligible in this range. It is a remarkable consequence of the high precision 
of our montecarlo estimates that we can still extract meaningful results by using only five values 

of A/3. 
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^3 extracted from ref. |34|]. Finally we can directly estimate the ratio by using 
the unbiased data for £ obtained in ref. [[M| and reported in tab. 9. The only 
problem is that these data are slightly less precise and that the lowest value 
of /3 is missing. The resulting estimate for the ratio: f = 1.029(11) is thus 
affected by a larger error. Also in this case we find a good agreement within 
the errors with our final result f ^ = 1.017(7) 


Table 11: Results of the fits according to eq. ( [Tip and (fP|) (see also the comment 
(a5) above). 



r x 

r € 

u 

r c 

£2 mi 

~X.red 

C.L. 

ao 

a\ 

0.52 

59% (2) 

3.5(8) 

46.9(5.4) 

0.43 

65% 

-0.04(50) 
-2.9(3.4) 

0.52 

67% 
4.8(1.2) 

0.01 

98% 

3.0(1.5) 
54(10) 

0.20 

90% 

0.24(8) 

final 

result 

4.75(3) 

1.95(2) 

14.3(1) 

3.05(5) 

1.009(5) 


5.1 Comparison with other existing estimates 


In tab. 14 we have compared our results with those obtained with other methods. 
Let us briefly comment on this comparison: 


bl] There are three possible approaches to the evaluation of the amplitude ratios. 
Montecarlo simulations (“MC” in tab. 14), low and high temperature series 
expansions (“HT,LT”, in tab. 14) and held theoretic methods. In this last 
case two different approaches are possible. The first one consists in looking at 
the e expansion of the £> 4 theory around four dimensions (“e-exp.” in tab.14). 
The second one consists in looking directly to the </> 4 theory in three dimen¬ 
sions (“d=3” in tab. 14). For a detailed discussion of these approaches see for 
instance ref. [35 . Let us also mention for completeness that an independent, 


interesting method to evaluate the ratio of specific heat amplitudes (which we 
do not study in this paper) by looking at the distribution of the zeroes of the 
partition function was proposed and applied to the Ising model in [RBI, R7I] . 


b2] The results of ref. [|. [|] for T x and Tg were obtained with a careful resummation 
of two loop e-expansions. On the contrary, the e-expansion for the exponential 
£ (which is needed to obtain the value of 2 nd which is reported in the last 

column of tab.14) is known only at one loop, hence the value f~/ f-, 2 nd ~ 1.005 
of [|J must only be considered as indicative. Later the e-expansion for T x was 
extended up to e 3 and the value of [fl] corresponds to the Pade resummation 
of such series. Recently in || this result has been further improved by using 
the parametric representation of the equation of state of the theory. 
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b3] The d — 3 approach originates from a suggestion of Parisi []T7|. While the 


results of || make use only of series expansions obtained in the symmetric 
phase of the theory, in [|7], || a three loop calculation, directly performed in 
the low temperature phase, was used. In this last case a crucial role is played 
by the low temperature renormalized coupling constant u evaluated at the 
critical point. We shall further comment on this point later. 

b4] The estimates for the amplitude ratios obtained with the use of low and high 
temperature series expansion reported in tab. 14 are mainly taken from ref ||. 
They were obtained with the use of IDA’s on the series reported in tab.12, 
where we used the standard notations: v = th(/3) and u = e -4/3 (not to be 
confused with the low temperature coupling constant in the d = 3 0 4 theory!) 
for high and low temperature series respectively. Recently these data have 
been reanalyzed in |L(| (using the same series) leading to essentially the same 
results. Notice however that in these last years new longer series have been 
constructed in the low temperature phase. These are the series that we have 
used in the previous sections to test our montecarlo results. It would be very 
interesting to see if these new series can lead to improved estimates of the 
amplitude ratios. In particular it would possible now to analyze also the ratio 
£/£ 2 nd which was previously unaccessible, since the two series for £ and for 
C, 2 nd start to be different only at the order u e 


Table 12: Some informations on the series used in ref.[9] 


ref. 

year 

observable 

length 


3S| 


1979 

HT/x 

v ls 


39 


1969 

HT/ £ 2nd 

v 12 


41 


1979 

LT/x 

u 20 


4C 


1975 

LT/ £ 2 nd 

u 15 


4C 


1975 

LT/£ 

u 7 


41 


1979 

LT/ m 

u 21 


Table 13: Some informations on the low temperature series used in this paper. 


ref. 

year 

observable 

length 


28 


1993 

LT/x 

u 32 


26 


1995 

LT/ £ 2n rf 

if 23 


27 


1995 

LT/ £ 

u 15 


28 


1993 

LT / m 

u 32 


b5] Some of the data reported in tab. 14 (those which are underlined) have been 
obtained combining separate amplitudes reported by the authors, thus their 
errors are most probably overestimated. 
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b6] It has been recently reported in ||43|| an estimate for the ratio obtained 


with a strong coupling expansion to 15th order of the correlation function 
G(t). 


The result is j- = 1.00023(5) which agrees with our bounds 1. < 


f -±- < 1.0006. 


f-\-,2nd 


Table 14: Results for the amplitude ratios reported in the literature 


ref. 

year 

method 

b'+ 

C- 

/+,2 nd 

u* 

c+ 

f- 

J — ,2nd 

n. 2 nd B2 

f— ,2nd 

ra 

1974 

e-exp 

~ 4.8 

~ 1.91 



~ 1.005 


i 


1985 

e-exp 

~ 4.9 






i 


1996 

e-exp 

4.70(10) 






a 


1987 

d = 3 

4.77(30) 



3.02(8) 



§ 


1996 

d = 3 

4.82(10) 






i 


1996 

d = 3 

4.72(17) 

2.013(28) 

15.1(1.3) 




i 


1989 

HT,LT 

4.95(15) 

1.96(1) 

14.8(1.0) 

3.09(8) 



42 


1993 

HT,LT 



14.73(14) 
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1994 

MC 

5.18(33) 

2.06(1) 

17.1(1.9) 

3.36(23) 


This work 

1996 

MC 

4.75(3) 

1.95(2) 

14.3(1) 

3.05(5) 

1.017(7) 


5.2 Comparison with an effective potential model. 


It has been recently proposed to study the critical properties of the three dimen¬ 
sional Ising model by constructing the effective potential of the corresponding quan¬ 
tum held theory. This effective potential is constructed simulating the model for 


various values of the external magnetic held. This program was carried on in [44 


for the high temperature phase of the model and was recently extended to the bro¬ 
ken symmetric phase in |4^|. The main result is that in the effective potential, 
besides the expected 0 4 term also a 0 6 term is present. It is interesting to test this 
model with our high precision results. Fortunately one of the values of f3 studied 
in |45|]: f3 = 0.2260 exactly coincide with one of our values thus allowing a detailed 
comparison. This comparison is reported in tab. 15, where in the third column we 
have reported the values of the various observables directly measured at (3 — 0.2260 
(tab.l of (|5|]) while in the last column we have reported the same observables ob¬ 
tained from what was considered in [|4l| as the most successful fitting procedure 
(tab.2 of [[45]). Notice that x is the inverse of V" and that u = 3 G in the notations 
of O ur results are in general one order of magnitude more precise than those 
151 


of ref. 


It is interesting to see that both m and x are i 11 rather good agreement 


with our data. The only strong disagreement is in the value of £2 n d and, as a con¬ 
sequence of this, in u. Most likely this disagreement is only due to the too small 
lattices studied in ref. (45] (the lattice size is reported in the last line of tab. 15) and 


does not imply that the approach proposed in j45| is wrong. 
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Tabic 15: Comparison with Tsypin’s results. 


Observable 

This work 

Tab.l of [0 

Tab.2 of 0 

m 

0.449984(16) 

0.44975(17) 

0.44975(17) 

X 

27.596(11) 

27.397(75) 

27.586(198) 

s2 nd 

3.0340(32) 

2.946(6) 

2.956(10) 

u 

14.64(5) 

15.90(9) 

15.84(28) 

L 

80 

30 

30 


5.3 Comparison with experimental data 

The experimental data reported in tab. 16 refer to the three most important ex¬ 
perimental realization of the Ising universality class, namely binary mixtures (bm), 
liquid-vapour transitions (lvt) and uniaxial antiferromagnetic systems (af). It is 
important to notice that these realizations are not on the same ground. Antifer¬ 
romagnetic systems are particularly apt to measure the C + /C- and f +, 2 nd/ f-, 2 nd 
ratios, while for the liquid vapor transitions the T c = R c /R | combination is more 
easily accessible. Finally, in the case of binary mixture all the three ratios can be 
rather easily evaluated. Even if obtained with very different experimental setups 
all these estimates qualitatively agree among them and this is certainly one of the 
most remarkable experimental evidences of universality. When looking in more de¬ 
tail at the various results one can see a residual small spread among them (even 
if in general the various estimates are compatible within the quoted experimental 
uncertainties). This spread is mainly due to the presence of correction to scaling 
terms whose amplitudes vary as the experimental realizations are changed and that 
are difficult to control. Thus some care is needed to compare these experimental 
data with theoretical estimates. The common attitude is to assume that the above 
systematic errors are randomly distributed and to take the weighted mean of the 
various experimental results. 


Table 16: Experimental estimates for some amplitude ratios 


exp. setup 

c± 

C- 

/+,2 nd 

c+ 

f — ,2nd 

fl.2n d B ' 2 

(bm) 

4.4(4) 

1.93(7) 

3.01(50) 

(lvt) 

4.9(2) 


2.83(31) 

(af) 

5.1(6) 

1.92(15) 


(all of them) 

4.86(46) 

1.93(12) 

2.93(41) 

This work 

4.75(3) 

1.95(2) 

3.05(5) 


Following this line we have reported in tab. 16 the weighted means (together with, 
in parenthesis, the standard deviations), of the experimental results reported in 0. 
In the first three rows we have studied separately the three different realizations of 
the universality class while in the fourth row all the experimental data at disposal 
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are analyzed together. In the last row we have reported our results. 

In the case of the f+^ n d/f-, 2 nd ratio ( for which, as we have seen, some of the 
present theoretical or montecarlo estimates disagree) we have listed, for a more 
detailed comparison, all the available experimental data in tab.17. In this table 
we denote with “N-H” the nitrobenzene - n-hexane binary mixture, and with “I- 
W” the one obtained by mixing isobutyric acid and water. A much more detailed 
account of the various experimental estimates can be found in |l;). 


Table 17: Experimental estimates for the { +,2nd ratio 

J— ,2nd 


Ref. 

year 

exp. setup 

f-\-,2nd 

J — ,2nd 


46 


1971 

(af),MnF 2 

1.7(3) 


47 


1972 

(af),FeF 2 

2.06(20) 


48 


1980 

(af),CoF 2 

1.93(10) 


49 


1983 

(bin),N-H 

1.9(2) 


56 


1986 

(bm),I-W 

2.0(4) 

This work 

1996 

MC 

1.95(2) 


6 Conclusions 


We have estimated various universal amplitude ratios in the case of the three di¬ 
mensional Ising model. Our final results are: 


A = 4.75(3) 


/+. 


2nd 


3 <7_ 


u 


■ 

= 14.3(1) 


= 1.95(2) 


/- 


C 4 


f 3 R2 v-' f 

J—,2nd J ~ > J+,2nd 


B 2 


f—,2nd 

= 3.05(5) 


= 1.017(7) 


(45) 


(46) 


Our results are in general in good agreement with other estimates of the same 
quantities obtained with held theoretical methods or with high/low temperature 
series. The main discrepancy that we have found is with the Montecarlo results of 
ref. m and with some of the results of ref. 1051. It must also be noticed that our 


result = 1.95(2) only marginally agrees with that of [§]: 2.013(28). However, 
as we mentioned above, this result depends on the value of the coupling constant u* 
which is an external input in the calculations of || . By plugging our value of u* in 
the perturbative expansion of [|J we find a lowering of the ratio with respect to j8|. 
This lower result: j + ’^ = 1.99(2) agrees not only with the strong/weak coupling 
result ||, but also with our estimate. Finally it is important to notice that our 
results are also in reasonable agreement with the experimental ones. 
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